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ABSTRACT 


The  problem  of  the  stress-induced  ordering  of  solute  atoms, 
the  so-called  Zener  effect,  is  considered.  In  this  work  the  effect 
of  the  difference  in  solute -solvent  atom  sizes  is  considered. 

Starting  with  this  size  factor  the  tensor  anelastic  strain  is  computed. 
Then  by  the  thermodynamic  theory  for  relaxation  processes,  due 
Meixner  and  others,  the  relaxation  strength,  is  calculated  as  a 
function  of  orientation.  The  results  of  the  calculation  suggest 
that  the  size  factor  is  not  necessarily  the  most  important  contri¬ 
bution  in  determining  the  relaxation  strength.  It  is  suggested 
that  electronic  effects  may  be  more  significant. 
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On  the  Zener  Effect  In  Ailidye 


I.  The  Size  Factor 

Introduction . 

Ordinarily,  an  externally  induced  change  in  state  of  a 
crystalline  solid  such  as  an  elastic  deformation,  magnetization  or 
electric  polarization  is  reversible.  However,  when  such  reversible 
processes  are  coupled  with  diffusion-type  mechanisms  the  atoms  of 
the  crystal,  assisted  by  thermal  fluctuations,  will  exceed  activa- 
tian  energy  barriers  and  migrate  to  neighboring  lattice  sites.  The 
previously  reversible  phenomena  are  now  retarded  and  damped;  for 
example,  the  strain  lags  behind  the  stress  field,  etc.,  vibrational 
energy  of  an  oscillating  system  is  transformed  irreversibly  to 
thermal  energy  at  a  finite  rate.  Such  phenomena  are  called  relaxa¬ 
tions.  Systems  which  show  a  specific  capacity  for  dissipating 
energy  internally  are  said  by  Zener  to  exhibit  internal  friction^^\ 

Now  there  are  a  number  of  micromechanisms  which  can  cmitribute 
to  the  internal  friction  of  a  crystalline  solid.  The  one  we  shall 
consider  in  this  paper  is  that  of  stress -induced  ordering  of  solute 
atoms  in  solid  solutions.  In  particular,  we  shall  consider  the 
effect  of  the  so-called  Size  Factor  of  the  solid  solution. 

If  we  assume  that  we  have  a  system  undergoing  periodic 
mechanical  vibrations  (e.g.,  a  torsional  pendulum),  which  were 
initiated  with  a  certain  amplitude  and  allowed  to  decay  freely, 
the  amplitudes  of  the  oscillations  will  die  out  exponentially  in 
time.  Ihe  internal  friction  is  then  usually  expressed  quantitatively 
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in  terms  of  the  logarithmic  decrement  6  of  the  oscillations  which 
we  define  by 

<f  -  (1) 

If  instead  the  system  is  driven  at  constant  an^litude  then 

(f  •  ^(AE/E)  (2) 

where  ^  is  the  energy  lost  to  heat  per  cycle  and  E  is  the 
maximum  amount  of  energy  stored  in  the  specimen  during  one  cycle 
of  oscillation. 

In  1948  Zener  proposed  the  following  dynamical  equation  in 
giving  a  general  phenomonological  description  of  anelastic  behavior, 

(T  S’  •  M^Ce-hTj:) 

Tf  ■  the  relaxation  time  of  stress  at  constant  strain, 

■  the  relaxation  time  of  strain  at  constant  stress,  and 

■  the  relaxed  modulus 
r 

For  those  processes  characterized  by  a  single  relaxation  time 
T  Zener  obtained  the  result 


/ 

TT 


tanf^» 


1  -  C*T)^ 


X  "  (T^  ^  phase  angle  by  which  the 
strain  lags  the  stress,  and  =  the  unrelaxed  modulus.  Q"^  is 
called  by  Zener  the  internal  friction.  If  we  plot  Q“^  versus  the 
frequency,  m  ,  we  obtain  a  bell-shaped  curve  whose  maximum  occurs 


-  3  - 

at  a  frequency  given  \uT  *  1,  see  Figure  1. 

In  general,  the  internal  friction  has  been  observed  to  be 
qviite  frequency  dependent,  usually  showing  successive  regions  of 
high  and  low  damping.  The  peaks  in  the  damping  behavior,  i.e,  of 
high  energy  loss,  are  called  relaxation  peaks  and  are  the  resultant 
effect  of  the  coupling  of  some  irreversible  process  in  the  material 
with  the  applied  stress.  In  our  case,  that  of  the  Zener  effect, 
the  internal  friction  peak  arises  from  the  coupling  of  the  stress 
with  the  diffusional  motions  of  atoms  in  substitutional  alloys  as 
they  become  preferrentially  ordered  in  the  stressed  crystal. 

In  situations  where  the  relaxation  processes  are  thermally 
activated  v/e  have 

T  *  t',exp(i^/Rr);  Ztfi  •  the  activation  energy  of  the 

relaxation  (5) 

so  that  we  may  conveniently  vary  the  relaxation  time  instead  of 
the  frequency  in  order  to  observe  a  peak.  Figure  2  shows  a  typical 
Zener  relaxation  peak  as  a  function  of  temperature  obtained  by 
Nowick  for  a  Ag-30/lo  Zn  alloy. 

There  are  two  parameters  which  characterize  a  relaxation  peak 
i)  the  relaxation  strength 
ii)  the  relaxation  time  T 

^  is  a  measure  of  the  peak  height  and  is  a  thermodynamic  quantity, 
depending  only  upon  the  initial  and  final  states  of  the  crystal. 
Following  Zener's  treatment  we  define  by 


(6) 
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T  is  a  kinetic  quantity  and  cannot  be  treated  by  equilibrium 
methods;  in  our  range  of  interest  ^  ■  the  activation  energies 
for  diffusion  in  the  alloy. 

Zener  has  described  the  internal  friction  peak  bearing  his 
name  as  follows:  Under  the  action  of  an  externally  applied  stress 
the  internal  stresses  tend  to  be  relieved  if  the  distribution  of 
atoms  (local  order)  is  altered,  thus  the  term  stress -induced 
ordering.  In  particular,  the  relaxation  is  due  to  the  reorienta¬ 
tion  of  pairs  of  adjacent  solute  atoms  into  preferrential  positions 
under  the  influence  of  an  applied  shear  stress.  Such  a  relaxation 
is  possible  because  the  size  difference  between  solute  and  solvent 
atoms  results  in  a  tetragonal  distortion  of  the  lattice  when  like 
pairs  of  atoms  are  similarly  oriented.  The  number  of  solute  atom 
pairs  was  assumed  by  Zener  to  remain  constant  throughout  the 
process.  This  model  seems  to  be  in  general  agreement  with  the 
result  that  relaxations  are  difficult  to  observe  in  systems  which 
show  a  small  size  factor. 

However,  as  pointed  out  by  LeClaire  and  Lomer, '  '  and 

Nowick  and  his  co-workers the  problem  is  of  considerably 
greater  complexity  than  indicated  by  the  simple  pair  model.  The 
configurational  arrangement  of  atoms  in  the  crystal  is  much  more 
complex  at  high  concentrations  (15-50%)  than  can  be  described  by 
a  simple  distribution  of  solute  atom  pairs.  No  account  is  taken 
of  the  possibility  of  short  range  order  or  clustering  being  present 
in  the  alloy.  Higher  order  neighbor  interactions  are  ignored  and 
in  general  the  number  of  pairs  of  solute  atoms  is  not  ccmstant. 


— I 
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Now  it  is  Itnown  that  alloys  while  in  the  stress -free  state 

shov;  isotropic  short-range  order  parameters.  On  the  other  hand, 

it  can  be  shown  tloat  the  local  order  in  a  stressed  solid  is 

directional.  In  fact,  the  short-range  order  is  a  symmetric  second 

rank  tensor,  te  are  then  required  to  specify  a  number  of  directional 

ordering  parameters  which  are  consistent  with  the  crystal 

neighbor  relations.  These  may  then  be  used  to  form  a  short- 

range  order  tensor  which  is  related  to  the  anelastic  strain.  To 

do  this  we  shall  follow  the  model  proposed  by  Chikazumi^  ^  ^  and 
(o') 

others^  '  '  in  treating  the  magnetic  annealing  behavior  of  ferro¬ 
magnetic  alloys.  In  making  this  calculation  only  the  atomic  size 
difference  is  taken  into  account.  After  doing  this  we  will  then 
calculate  ^  and  its  behavior  as  a  function  of  orientation, 
ten^rature,  etc. 

Given  a  crystal  of  A  and  B  type  atoms  having  N  sites 
there  are  nearest  neighbor  pairs  of  all  types.  We  write  this 
as 

[aa]  +  [ab]  +  [ba]  f  [BB]  *  ^  -  nN  (7) 

where  ■  the  number  of  AA  pairs,  etc.;  n  ■  JjZN  ■  the 

number  of  independent  orientations  possible  for  each  pair. 

Considering  the  total  number  of  nearest  neighbor  pairs  aligned 
the  direction  we  have 

[aa]  .  +  M  +  Cmlj  +  Cb]3  j  «  M  (8) 


(9) 
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If  one  furtiier  defines  a  concentration  parameter  C  to  be 

C  -  (  {a]  -  [f]  )/N 

and  a  directional  order  parameter,  for  nearest  neighbor  pairs 

along  the  j  direction  to  be 

-  r{  [“]  J  ♦  W  j  -  M  3  -  [»?  j} 


one  obtains 


and 


M  j  .  If  (1  t  2C  -or^) 
[bBJ  j  .  I  (1  -  2C  -  01^) 


From  this  result  one  can  then  write  down  the  average  atom  spacing 
along  the  j  direction  of  the  stressed  crystal  as 

dj  =  jr  [(r.  ♦  Ajj)  [m]  .  ♦  (r.  *  AjjX  ^  ♦  [»] 

-  (r.  *  o^rgg)  ^'j 

where  r*  ■  the  average  atom  diameter  determined  from  the  lattice 
parameter.  By  substitution  tdiis  result  reduces  to 


dj  *  r.  -  V(l+2C-aj)<fr^  4  %(l-a^)cfr^g  f  ^^(1-20-0^  )<yrgg 

The  average  spacing  d^  along  the  j  direction  of  the  stress - 
free  crystal  is  given  simply  by  substituting  the  isotrpplo  value. 


a,  for  the  short  range  order  parameter  a. .  The  linear  dilatation 
along  the  direction  due  to  the  stress 'induced  ordering  of  atoms 
is  then  given  by 

(d^-3j)/c[.  -  -e(Oj-a)/2r.  (12) 

ft  sm 

\7here  9  s  related  to  deviations 

from  Vegard's  law  for  lattice  parameters.  When  9*0,  Vegard's 
rule  is  obeyed.  To  first  order  we  may  also  assume  that  the  linear 
dilatations  along  directions  normal  to  the  j  direction  due  to 
ordering  along  direction  j  are  zero. 

We  now  write  the  total  contribution  to  the  strain  due 
to  ordering  along  the  j  direction  as 

•21  dU^  A  A 

C.  z  —  0*  cr  (13) 

J  do 

-0  — 

where  du^/dc  *  (a^-  at)  and  is  a  cartesian  coordinate 

along  direction  j;  ^  is  a  unit  vector  along  the  j  direction 
which  we  may  write  as  a  r  ^jy^  ^  ^jz^'  ^ji 

the  direction  cosines  of  o'  relative  to  the  i^*'  crystal  axis. 

Taking  account  of  the  contributions  to  the  strain  due  to 
ordering  along  all  of  the  j  directions  we  have 


e 


rHfJ 


j  -  1 


(14) 


is  a  short'Zange  order  tensor  \^ch  we  define  by 


8 


yik 


(cc.  -  o) 


(IS) 


In  the  case  of  the  body-centered  cubic  and  face-centered  cubic 
structures  theA. a  £(3)”^,  and  0  and  ^(2)“^  respectively.  In 
terms  of  this  result  we  should  note  that  when  Vegard’s  Law  holds 
(0  ■  0)  there  will  be  no  anelastic  strain  associated  with  stress- 
induced  changes  in  order  which  arise  from  the  size  factor  of  the 
alloy.  This  is  in  agreement  with  the  conclusions  of  LeClaire  and 
Lomer^  ^  ^  based  on  their  earlier  theory  of  directional  ordering 
effects  in  internal  friction. 


Calculation  of 

Now  to  calculate  ^  we  must  determine  the  Helmholtz  free 
energy  F  as  a  function  of  (C,  0,  T,  o^.).  Then,  by  the 
theiTOodynounic  theory  of  relaxation  processes,  we  calculate  ^ 
as  a  function  of  orientation,  temperature,  etc. 

Following  OeDonder,  Meixner  and  others^® ^  we  have  for  the 


case  of  isothexmal  relaxations 
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Z  §3 

j 


V  L 


<r>2^ 

j 


(16) 


where  both  a.  and  A.  are  linear  functions  of  (T,  C  and  T.  A.  is 
J  J  3 

the  Chemical  Affinity  corresponding  to  o(^.  and  is  defined  by 


-A.  s  dF/do. ;  A.  >  0  at  equilibrium. 

*0  J 

Dividing  each  j-term  on  both  sides  of  the  equation  by  <  R  > 
(dAyd<y)g^j  we  obtain  on  re-arranging  terms 

I  •  «r  [  II  (17) 

j  ^  j  J 

Comparing  this  result  witii  Zener’s  equation  (3)  we  may 
write  that 

and  _  . 

^  =  Knr^ 

And  from  this 


3 

3 


<  R  >  represent^  an  appropriate  Arxheniua  rate  factor. 


since  M  ■ 

u  r 


From  the  thermo<tynamlc  relation 


we  may  then  show  that 


^  ‘ 


Cd»3/a<»j)e,T  - 


06) 


a?) 


vftiich  is  a  result  similar  to  that  of  Zener. 

In  general  we  can  break  the  Free  Energy  up  into  three  terms 

(1)  Internal  energy  including  chemical  interactions 
and  electronic  effects,  U. 

(2)  Strain  energy  contributions,  W. 

(3)  Configurational  entropy  terms, 

Following  Iwata^^  ^  vve  consider  that  the  Qoasi-Chemical  Approximation 
applies  and  vie  then  write 


j 


where  V  =  ^  ”  '^AB'  >  0  then  we  tend  to  get 

ordering;  if  <^<  0  we  tend  to  get  clustering. 

Idnear  elasticity  theory  gives  the  strain  energy  to  be 


u 


w  =  ^  ‘^ik^ik  •  (19) 

i,k  i,k 

C  7  ^ 

Iwata'*^  *'  has  also  calctilated  the  configurational  entxopy 
for  nearest  neighbors  by  the  self-consistent  me^od  of  lakagi^^  ^ 
which  is  obtained  from 


■  !<•  log^w 
conf  ^e 

where  W  *  (nN)!  PjPjj/  [^1^.1  [AB]j![  BAj^!  [bb]^  !  . 

the  thermodynamic  probability. 

for  the  case  of  a  body-centered  cubic  lattice.  Pj  and  are 

the  probabilities  for  all  consistent  occupations  of  the  sub-lattices 
I  and  II.  A  similar  result  is  obtained  for  the  face-centered  cubic 
lattice,  which  has  four  sub-lattices.  (The  equilibrium  values  of 
the  a.  may  be  obtained  from  the  condition  A.  ■  0).  The  result 
tlien  gives  the  configurational  entropy  as 

Sconf  *  where  is  a  long  range  order  parameter. 

If  we  consider  the  case  v/here  ;•/  •  0,  and  with  the  Free 
Energy  we  have  calculated,  F((5,e,T,  x^,a^)  ■  U  -  W  -  we 

may  then  evaluate  A.  and  hence  We  obtain  the  resvilt  tiiat 


-A. 


aF 


kTN 


Q(a.,x„) 


i,k 


3 


(20) 


12 


Q(^j>  ^  function  of  the  directional  ordering  parameters 

a.,  and  the  atom  fractions  of  the  constituents,  x^.  B/  substitution 
into  equation  (17)  we  may  then  calculate  the  relaxation  strengtdis 
of  BCC  alloys  associated  with  the  various  components  of  the 
anelastic  strain.  We  give  tdve  results  for  tensile  and  shear 
relaxations  along  the  <  100  >  and  <  111  >  directions  respectively. 


^,100 


T  -  Ej^ju|8*/9Nk?rJ  ’ 


(21) 


^,111 


i'  (e/2r,)^ 

/t  -  4Ett,qV4NI^? 


3(0/18r.)2 

T-4ET,,0^NlJ"18r 


are  the  relaxation  strengths  along  the  <  100  >  and  <  111  > 
directions  associated  with  the  Young’s  Modulus. 

is  a  function  of  the  ordering  parameters,  o^.,  and  the  concentra¬ 
tions,  x^.  Ej^qq  and  E,.,^  are  the  unrelaxed  Young’s  moduli  for  the 
<  100  >  and  <  111  >  directions  respectively. 

4Gj^009='/llto’9r; 

'  T  -  2Gjjpe‘'/l(W’9r? 

and  _  (22) 

4G^j^O^/Nk(^7r^ 

*  T  - 

are  the  relaxation  strengths  in  the  <  100  >  and  <  111  > 
directions  associated  with  the  shear  modulus  G.  G^qq  and  G.,^., 
are  the  unrelaxed  shear  moduli  for  the  <  100  >  and  <  111  > 
directions  respectively.  A  similar  set  of  relations  may  be  obtained 


1 
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for  the  analogous  relaxation  strengths  in  FCC  alloys. 

Discussion . 

If  we  consider  the  results  of  the  calculation  certain  conclusions 
may  be  drawn.  (1)  equation  (17)  gives  a  general  expression  for  an 
anisotropic  ^  which  is  independent  of  a  specific  model  for  the 
relaxation.  (2)  the  Size  Factor  contributions  to  the  relaxation, 
given  by  (21)  and  (22),  predict  results  which  contradict  the 
directional  dependence  tliat  is  observed  experimentally;  see  Figures 
3  and  4.  (3)  a  relaxation  of  a  pure  hydrostatic  pressure  is 

predicted  by  equations  (17)  and  (19)  This  agrees  with  the  predicticms 
of  LeClaire  and  Lomer^  (4)  the  tempezeture  dependence  predicted 
by  equations  (17)  and  (20)  for  is  that  which  was  first  predicted 
by  Zener^^^  and  observed  experimentally  by  Li  and  Nowick^^\ 

The  anisotropy  ratios  ^,ioo/^,lll  ^,lll/^,100» 

calcxilated  from  equations  (21)  and  (22),  give  the  wrong  directional 
dependence  when  corapared  with  experimental  results  for  A 

similar  calculation  made  for  FCC  systems  is  likewise  in  error  in 
the  same  sense.  The  anisotropy  in  ^  is  the  same  for  FCC  as  for 
BCC  alloys^^^^^^  Our  conclusion  is,  therefore,  that  the  atomic 
size  difference  cannot  account  for  these  relaxations.  Furthermore, 
if  vie  substitute  the  appropriate  values  of  E,G,  0,  jp*  and  r, ,  into 

equations  (21)  and  (22)  the  calculated  values  obtained  for  ^  and 
-4 

^  are  ■  10  and  are  at  least  an  order  of  magnitude  smaller 
than  the  observed  relaxation  strengths  given  in  Figuzes  3  and  4. 

The  con^sition  dependence  of  enters  thxougtt  tite  dimensionless 
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number  ^  ,  which  is  approximately  proportional  to  the  sctuare  of 
the  solute  concentration  in  the  dilute  case. 

The  lack  of  agreement  between  experiment  and  predictions  based 
on  nearest  neighbor  pair-reorientation  theories  was  explained  by 
Nowick  et  as  being  due  to  neglect  of  next  neighbor 

pair  terms.  To  check  this  possibility  our  calculation  was  modified 
to  consider  next  neighbor  ordering  directions  also.  Our  results 
did  not  show  any  change  in  the  character  of  the  directional 
dependence  and  hence  no  improvment  in  agreement  with  experiment.  It 
does  not  appear  that  any  further  generalization  of  the  atomic 
size  factor  contribution  will  lead  to  any  improvement.  Attempts  to 
include  chemical  interactions  v;ith  the  size  effect  contribution 
have  also  met  with  little  success'  In  our  view,  the  reason 
for  this  is  that  the  relaxation  is  closely  related  with  perturbations 
of  the  electronic  structure  of  the  alloy.  This  is  because  of  the 
known  dependence  of  the  electronic  structure  on  the  degree  of 
short-range  order  in  the  alloy.  Tne  most  likely  contributions 
are  those  involving  interactions  between  the  Fermi  surface  of  the 
alloys  and  their  average  Brillouin  zone  boundaries.  We  are 
investigating  this  problem  at  the  present  time  and  will  discu.ss 
it  in  a  later  report. 

The  author  would  like  to  express  his  thanks  to  the  ONR  for 
their  generous  support  of  this  work. 
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Figure  Captions 


Figure  1.  Schematic  representation  of  a  relaxation  peak  In 
the  internal  friction.  The  S>shaped  curve  gives  the 
frequency  dependence  of  the  associated  elastic  modulus. 

Figure  2.  The  internal  friction  peak  for  two  Ag-Zn  alloys  of  30.25% 
and  16.5%  Zn  composition  according  to  NowiCk  et  al.(4) 

Figure  3.  The  relaxation  strengths  ^  and  for  FCC  single  crystal 
wires  (4)  of  different  orientations.  The  alloys  were 
Ag-26  at  .%  Zn  by  composition.  The  function  n  defines 
the  orientation  of  the  wire  axis  where 

>^2*  direction  cosines  of  the  axial 

direction  relative  to  the  <100>  directions  of  the 
crystals.  A  similar  behavior  is  observed  for  BCC  Li -Mg 
alloys . 

Figure  4.  Orientation  dependence  of  ^  for  Ag-Zn  alloys  (4) 
of  various  conqwsitions . 
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